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Abstract 

We have calculated all contributions to the energy levels of parapositronium at order ma 1 coming from vacuum polarization cor¬ 
rections to processes involving virtual annihilation to two photons. This work is motivated by ongoing efforts to improve the 
experimental determination of the positronium ground-state hyperfine splitting. 
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1. Introduction 

Positronium, the electron-positron bound state, is a particu¬ 
larly simple and interesting system. The constituents of positro¬ 
nium have no known internal structure. The properties of 
positronium are governed almost completely by QED-strong 
and weak interaction corrections are below the level of current 
interest due to the small value of the electron mass. On the 
other hand, some features of positronium tend to complicate 
the analysis compared to, say, hydrogen or muonium. The no¬ 
recoil approximation is not relevant for positronium-the mass 
ratio for positronium takes its maximum value of one. Also, be¬ 
cause positronium is composed of a particle and its antiparticle, 
it exhibits real and virtual annihilation into photons. The states 
of positronium can be taken to be eigenstates of the discrete 
symmetries charge parity and spatial parity, making positro¬ 
nium useful in searches for new, symmetry breaking interac¬ 
tions. Because of its unique properties and accessibility to high- 
precision experiments, positronium is an ideal system for tests 
of the bound state formalism in quantum field theory and for 
searches for new physics in the leptonic sector. 

Since its discovery in 1951 01, positronium has been the ob¬ 
ject of increasingly precise measurements of the ground state 
hyperfine splitting (hfs), orthopositronium (spin-triplet) and 
parapositronium (spin-singlet) decay rates and branching ratios, 
n = 2 fine structure, and the 2 S - IS interval. This progress is 
reviewed in Refs QS&lSH&SGl with citations to the 
original literature. The most precise hfs measurements were 
performed by two groups in the 70s and early 80s [11, 12 , 13]: 

AEfBrandeis) = 203 381.5(1.6)MHz, 

AE(Yale) = 203 389.10(74 )MHz. (1) 

Both of these results are based on the observation of Zeeman 
mixing of ortho and para states in the presence of a static mag¬ 
netic field. More recently, a great deal of work has been done 


both with the Zeeman approach and with other indirect and di¬ 
rect methods of measurement mmEESEiiiiiii]. a 
new high-precision measurement utilizing the Zeeman method 
with improved control of systematics was recently reported 1 22tl 


AE(Tokyo) = 203 394.2(1.6) stat (1.3) sys MHz. 


( 2 ) 


Theoretical work on the positronium hfs involves calculating 
the energy splitting by use of bound-state methods in QED. The 
principal modern approach involves the definition of an effec¬ 
tive non-relativistic theory through matching with full QED fol¬ 
lowed by a bound-state perturbation calculation in the effective 
theory. The result has the form of a perturbation series in the 
fine structure constant a augmented by powers of C — ln(l/a). 
This series has the form 


A E = ma 


{Co + C\ — + C2\a 2 £ + C 20 j 

(;)’*■•) <3, 


a 3 0 a 3 

+C 32 { + C31 { + C30 

n n 
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where m is the electron mass. The coefficients Co - C31 are 
known analytically as reviewed in Ref. |23!]. The most recent 
result, 

17 217 

C 31 =- y ln2 + —, (4) 

was obtained by three groups in 2000 |0 25, 2f|. The numer¬ 
ical value of the theoretical prediction, including terms through 
C31, is 

AE(th) = 203 391,69MHz (5) 

with an uncertainty due to uncalculated terms that has been esti¬ 
mated as 0.16 MHz 1 25 1, 0.41MHz [ 0 ], or 0.6 MHz ll23ll . The¬ 
ory and the older experiments are separated by 2.6cr and 3.5<x 
in terms of the experimental uncertainties, but theory and the 
new experiment are consistent with one another. 

The naive size of 0(ma 7 ) corrections is only ma 4 (a/ 7r) 3 = 
439kHz , but contributions as large as several tenths of a 
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MHz have been found coming from “ultxasoff ’ energy scales 
1 271 i28il . Additional Oinitr 7 ) contributions have recently been 
obtained [23, 2Sj 30[ 31]. The present work is a contribution to 
a systematic calculation of all corrections at 0 (ma 7 ) begun in 
anticipation of yet more precise measurements of the positron- 
ium hfs. 

Contributions to the positronium hfs at 0(ma 7 ) can be clas¬ 
sified as either annihilation or exchange depending on the pres¬ 
ence or not of virtual annihilation e + e~ —> ny —> e + e~ in 
the description of the process. Among annihilation contribu¬ 
tions, ones that involve virtual annihilation to an odd number of 
photons only affect orthopositronium according to charge con¬ 
jugation symmetry, while ones involving an even number of 
photons only affect parapositronium. We consider here two- 
photon-annihilation processes affecting parapositronium and 
focus specifically on processes containing a vacuum polariza¬ 
tion correction to one or both of the annihilation photons. This 
set of contributions forms a gauge invariant set, and further¬ 
more is insensitive to the particular bound-state formalism used 
in its evaluation, and consequently it comprises a reasonable set 
of contributions to be evaluated in isolation from other types 
of terms. Most of the two-, three-, and four-photon annihila¬ 
tion contributions have non-vanishing imaginary parts, as can 
be seen from Cutkosky analysis ll32ll . However, the vacuum po¬ 
larization function vanishes for small k 2 (where k 1 ' is the photon 
momentum), so the vacuum polarization corrections discussed 
here vanish for annihilation to on-shell virtual photons. These 
contributions to the energy shift A E are purely real and that fact 
simplifies their evaluation considerably. 


2. Pure vacuum polarization corrections 


The “pure vacuum polarization” corrections involving either 
the two-loop vacuum polarization function or a product of two 
one-loop functions are shown in Fig. [T] The effect of a vacuum 
polarization (VP) correction is to modify a photon propagator 
according to 


1 1 1 ^ 1 
1 , n t w ~ “ ~ ~ n R (/r)— + ■ ■ 

1 + n R (p 2 ) p- p- p- 


(6) 





Figure 1: The three types of 0{ma 7 ) pure vacuum polarization corrections in 
the two-photon-annihilation channel. Graph (a) represents the contribution of 
the one-photon-irreducible two-loop vacuum polarization function. This func¬ 
tion is represented by a single diagram, but there are two additional contribu¬ 
tions that can be described as self-energy corrections inside the vacuum po¬ 
larization loop that are not shown. Graph (b) gives the one-photon-reducible 
contribution. Graph (c) shows a final two-loop vacuum polarization correction, 
one with a one-loop correction on each virtual photon. Graphs with crossed 
photons are not displayed-the crossed photon graphs give energy contributions 
equal to those of the graphs that are displayed. In addition, the contributions of 
graphs (a) and (b) must be doubled since the correction could occur on either 
virtual photon. 


The energy shift due to the vacuum polarization graphs con¬ 
sidered in this article is insensitive to the particular bound state 
formalism we employ because the energy and momentum val¬ 
ues that contribute are purely “hard”-of order in, not ma or 
ma 2 . We choose to use the formalism of Ref. [37|], in which 
the energy shift is an expectation value 


A E = i'VSK'V. (9) 

The two-body states VP, 'F carry information about the positro¬ 
nium spin state X- 


where the renormalized scalar vacuum polarization function 
can be expressed in a spectral form as 


¥ 


0 X 
0 0 


T" 


00 


0 0 


or 


( 10 ) 


n r (p 2 ) 


-f 


dvg(v) 


p 2 — 4m 2 /(l — v 2 ) 


(7) 


At one-loop order the spectral function g(v) takes the form 


gi(v) = - 


v 2 ( 1 - v 2 /3) a 

1 - v 2 n 


( 8 ) 


The two-loop function has been given by Kallen and Sabry |33], 
Schwinger 13411 . and in a form using “standard” notation for the 
dilogarithm function ll35ll . by Eides, Grotch, and Shelyuto |36]. 
One must note that the Kallen-Sabry form is for the reducible 
vacuum polarization function (Figs, la plus lb), not the irre¬ 
ducible function of Fig. la alone. The Schwinger and Eides- 
Grotch-Shelyuto forms are for the irreducible function. 


where x i s the two-by-two two-particle singlet spin matrix 
X — 1 / V2. The positronium states, in this approximation, have 
vanishing relative momentum, and 0 0 = y/m 3 a 3 /( 87 m 3 ) is the 
wave function at spatial contact for a state of principal quantum 
number n and orbital angular momentum f: = 0. 

The explicit expression for the energy shift for the vacuum 
polarization diagram of Fig.QJ is 


A E 


1 a 


x tr 
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*00 J 


crp , „ , —i 

1 (~n«(/? 2 )) 


x tr 


^ 0, 

\i-iey v ) 


( 2 n ) 4 

i-ieY 1 ) 


p 2 (P - p ) 2 
- (~iey v ) 1 

— n 1 — m J 


y(P /2 - p) - m 


y(P /2 — p) - m 
(~ iey> J ) 


1° 

A] 

io 

ojJ 


] (ID 


2 





































where P = (2m, 0) is the positronium 4-momentum in the 
center-of-mass frame, p is the 4-momentum of the vacuum po¬ 
larization corrected photon, the initial (-1) is a fermonic minus 
sign, and the factor of 4 accounts for the graph with crossed 
photons and for the fact that the vacuum polarization correction 
could act on either photon. After evaluating the traces and some 
simplifications, the energy shift above takes the form 

AF = — f d4p 2p 2 n g (p 2 ) 
lo n J in 2 p 2 (p - P) 2 (p 2 - p ■ P) 2 

for the ground state (n =1). The integral over p can be eval¬ 
uated either using Feynman parameters or by the method of 
poles-closing the p° contour with a half-circle of infinite ra¬ 
dius in the upper or lower half plane. In either case, the result 
for the energy shift becomes 

4£ = I ,rf ' , d^t ,2 - v)<1+ ” )2l "( i ir) 

-(2 + v)(l-v) 2 ln(—^—) 

+2v 3 ln v + 2v(l - v 2 )}^-. (13) 

’ 7T 

For the one-loop vacuum polarization of ([Sj the energy shift is 


mar 1 

A Eyp-I = Iyp — , Iyp = ——£(2) (14) 

n~ 6 

in accord with the result of [38]. The two-loop irreducible vac¬ 
uum polarization contribution of Fig. la is 


AE la 


{——£(4) + —4(2) In 2 2 - — In 4 2 
l 2A h 3 18 

4 161 39 

— a A + -4T(3)-4(2) In 2 

3 96 8 


475 43 1 mop 

+ T92« 2) -96| — 


(15) 


where 4 is the Riemann zeta function and a A = Li4(l/2). We 
didn’t actually “do” the integral for A E\ a . Rather, we obtained 
a numerical result to high precision (100 digits) and used the 
PSLQ algorithm l39ll to obtain (HD . 

For the reducible contribution of Fig. lb we were able to ac¬ 
tually perform the exact integral. We used the spectral function 


gib(v) = 


2v 2 (1 - v 2 /3) i 8 v 


1 


{- 

l 9 3 


(1 



1 + v 
1 - v 



extracted from the work of Kallen and Sabry [S331 in (IT3l > to 
obtain 


A E xb = {—4(3) — -4(2)ln2+ — 4(2)- — (17) 

>6* 3 s 36 fe 108 ' tP 


Finally, for the term of Fig. lc with a one-loop vacuum po¬ 
larization correction on each annihilation photon, we were not 
able to achieve an exact result. The numerical value of this con¬ 
tribution is 


AEic = -0.045140511436— t-. (18) 






Figure 2: The four types of two-loop corrections to the two-photon-annihilation 
channel containing a one-loop vacuum polarization function. Graph (a) con¬ 
tains a self-energy correction. Graph (b) displays a correction having a vertex 
correction on the same photon that is modified by vacuum polarization (this is 
“Vertex Type A”). Graph (c) shows a correction where the vertex and vacuum 
polarization affect different virtual photons (this is “Vertex Type B”). Graph 
(d) contains a ladder correction. Each of these graphs must be multiplied by a 
factor of eight to account for equivalent graphs involving crossed annihilation 
photons, VP corrections on either photon, and the self-energy, vertex, or ladder 
corrections appearing at various positions in the diagram. 


3. Terms involving one-loop vacuum polarization 


Additional vacuum polarization corrections are shown in 
Fig. 2. They contain a one-loop vacuum polarization part com¬ 
bined with self-energy (SE), vertex, or ladder corrections to 
the bare two-photon-annihilation process. These contributions 
form a set that is best calculated together because, after renor¬ 
malization of the self-energy and vertex parts, each of these 
graphs contains an infrared divergence that only cancels when 
all four are summed. 

Our calculational method is to start with a form like (ED con¬ 
taining the one-loop vacuum polarization correction of ({7} and 
® and tack on a self-energy, vertex, or ladder correction. Con¬ 
venient forms for the one-loop self-energy and vertex correc¬ 
tions were given in |40[] ■ The renormalized self-energy function 
corrects the bare electron propagator according to 

—-^ -{^! + 5 2 (/z)}— i- (19) 

yp — m 1 yp - m 

where 



Slip) = J dxdufst—^ (20b) 
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with fsE = - 1 /( 2 u) and 


N(p) = {2m - (1 - x)yp] (yp + m) (21a) 

t t jwm 2 

Dip) = p -m 2 -- -—. (21b) 

(1 — x)u 

(All parametric integrals in the self-energy and vertex correc¬ 
tions run from 0 to 1.) A photon mass A was introduced as 
an infrared regulator in the course of on-shell renormalization. 
The vertex correction has the form 

/ -» + V^p',p) + V%(p',p)} (22) 


where p and p' are the incoming and outgoing electron mo- 


menta and 



v? 

= y(- in (-)-j), 

\ \ m ; 4 / 

(23a) 

V^ip'.p) 

= f dxdu ——, 

J AH 

(23b) 

V?ip',p) 

= / f Jxdudz-*"-™ 2 ). 

J 2 H 

(23c) 

The parametric functions are 


H = 

(1 - x)[u(m 2 - p' 2 ) + (1 - u){m 2 - p 2 )] 
—xu{ 1 — u)ip' - p) 2 + xm 2 , 

(24a) 

H = 

xm 2 + z(H - xm 2 ). 

(24b) 

W = 

y A (yip' + Q) + m) y (yip + Q) + m) y A 

(24c) 


where Q = -x(up' + (1 - u)p). 

The self-energy contribution is pictured in Fig. 2a. The di¬ 
agram shown must be multiplied by eight to account for the 
two places where the SE correction could occur, the two pho¬ 
tons that the VP could correct, and the two types of graph (with 
uncrossed and crossed photons). All eight contribute equally. 
There are two parts to the SE contribution: S \ and Slip) of 
( IT9t . The S i contribution is 2 S\Iypima 1 /n 3 ) where Iyp is the 
0(tna 6 ) VP correction of ( 1 1 4b , For the Slip) contribution we 
could preform the d 4 p integral as a whole using Feynman pa¬ 
rameters, or do dp 11 first by poles followed by d 3 p —> Anp 2 dp 
numerically. Both results are shown in Table Q] followed by 
their total. 

There are two classes of vertex corrections: ones with the 
VP and vertex corrections affecting the same photon (type A, 
Fig. 2b) or affecting different photons (type B, Fig. 2c). Each 
has three parts as specified in (l22l ) and < 1 23 b , The various vertex 
correction integrals were performed both by poles and parame¬ 
ters, with results recorded in Table [I] 

The ladder correction shown in Fig. 2d gives an energy shift 
that looks like ifTTb except that the right-hand trace of (ITU is 
replaced by 


tr[] -> 


x 

x 




(2 7 T ) 4 
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y(P/2 + q - p) - m 


y(—P/ 2 + q) — m 
(-*'<?/) 
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Table 1: Contributions to the p-Ps energy levels coming from one-loop vac¬ 
uum polarization (VP) corrections combined with self-energy (SE) (Fig. 2a), 
vertex (Figs. 2b and c), and ladder (Fig. 2d) corrections to the two-photon- 
annihilations graphs. For type A (type B) vertex corrections, the VP and vertex 
corrections act on the same photon (different photons). The separate contribu¬ 
tions to the self-energy, etc., parts are shown along with their totals, all in units 
of ma 1 171^. The infrared divergence is contained in L = In (A/m) Iyp. 


Term 

Poles Result 

Parameters Result 

SE: 1 

SE: 2 

SE 

2L~\C(2) 

-0.76281(16) 

2L- 1.03697(16) 

2L - A(2) 
-0.7629639(4) 

2L- 1.0371196(4) 

Vertex A: 1 

Vertex A: 2 

Vertex A: 3 
Vertex A 

- 2 L+&X 2 ) 

-0.2272831(5) 

0.1876296(12) 

-2L +0.6457357(13) 

- 2 L+ff( 2 ) 

-0.2272846(14) 

0.1876304(19) 

-2 L + 0.6457350(24) 

Vertex B: 1 

Vertex B: 2 

Vertex B: 3 

Vertex B 

-2L+±((2) 

0.0693357(4) 

0.2046073(11) 

-2L +0.9593322(12) 

- 2 L+^( 2 ) 

0.0693348(6) 

0.2046053(16) 

-2L +0.9593293(18) 

Ladder: 1 

Ladder: 2 

Ladder: 3 

Ladder 

2 L+K( 2 ) 

0.3124063(6) 

-0.5448084(11) 

2L +0.3159093(13) 

2 L+^( 2 ) 

0.3124057(21) 

-0.5448063(10) 

2L +0.3159108(24) 


where the factor of 2 comes because the ladder correction could 
occur on either side of the diagram. Again an infrared diver¬ 
gence, here arising from the binding singularity, is regulated 
by including a photon mass A. The contribution of (l25t can be 
written as 


- 2(47tq') 2 


/ 


d A q Niq) 
(2n) 4 D(q)Ziq) 


(26) 


where D{q) = ((P/2 - q) 2 - m 2 )((P/2 + q) 2 - m 2 )(q 2 - A 2 ), 
Z(q) = ((P/2 + q - p) 2 - m 2 ), and 


Niq) 


tr[/ iyi-P /2 + q) + m) y v iyiP /2 + q - p) + m) 
Xy^{y{P/2 + q) + m)y^ 0 *J], (27) 


We isolate the small-g singularity by use of the decomposition: 


Niq) 

D{q)Z{q) 


NiO) N{ 0)/ 1 1 \ 

Diq)Z( 0) + WrW) - w 
Njg)-N( 0) 

Diq)Ziq) 


(28) 


consisting of terms 1, 2, and 3, which we evaluate in turn. The 
first term in (l28l > contains the infrared singular part of the q 
integral, but otherwise is simply proportional to Iyp because 


WO) 

Z(0) 


, 2 tr[y v 

L yiP/2 - p) - m 





(29) 


is proportional to the right-hand trace of (1111) . The IR-singular 
binding integral is Hit 


f 


d 4 q — m 2 
in 2 Diq) 


In(-) - 1 + O(-). 

' m / ' m' 


(30) 
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Table 2: Positronium energy level corrections at ()(map) coming from vacuum 
polarization corrections to two-photon-annihilations graphs. The tabulated re¬ 
sults are contributions to I where A E = (map /7t 3 )I. The poles and parameters 
results from Table|T|were combined to give the results shown here. The infrared 
divergence contained in L = ln(/l/m) lyp cancels in the sum. 


Term 

Diagram 

Energy shift 

irreducible 2-loop VP 

la 

-0.9205630 

reducible 2-loop VP 

lb 

-0.0712481 

1-loop VP on each 

lc 

-0.0451405 

VP-SE 

2 a 

2 L- 1.0371196(4) 

VP - vertex A 

2 b 

-2L +0.6457355(12) 

VP - vertex B 

2 c 

-2L + 0.9593313(10) 

VP - ladder 

2 d 

2L +0.3159096(12) 

total 


-0.153095(3) 


The mn/A term represents the part of the ladder correction that 
comes from exchange of a Coulomb photon and so is just part 
of the binding that created the bound state in the first place, and 
so must not be included again here [42]. Or, to put it another 
way, the mn/A term is the piece that is removed in the matching 
calculation if an effective field theory formalism had been used 
[43, 41]. So term 1 contributes 2(ln (A/m) - 1 ) Iypima 1 /n 3 ) to 
the energy shift. Terms 2 and 3 of the ladder correction were 
constructed to be IR safe and can be evaluated directly using 
A —» 0. (It is useful to average over the directions of q when 
performing the d 3 q integration as described in Hi .) The re¬ 
sults for the various parts of the ladder correction are reported 
in Table 1. 


4. Results and discussion 

The results for all VP corrected two-photon-annihilation 
graphs at 0(ma 7 ) are shown in Table[2] We see that all infrared 
divergences cancel in the sum. The ground state parapositron- 
ium energy level correction due to vacuum polarization effects 
in the two-photon-annihilation channel is 

innA 

AE = -0.153095(3)—(31) 
n 3 

The numerical value of this correction is small due to large can¬ 
cellations among the various terms. This parapositronium shift 
corrects the hfs (spin-triplet minus spin-singlet) by the amount 

A£ hfs = 0.67 kHz. (32) 

Higher S states with principal quantum number n have correc¬ 
tions that are the same as (OTb except divided by the n 3 factor 
that comes from the square of the wave function at spatial con¬ 
tact. States with C 4 0 are not corrected by the terms consid¬ 
ered here at order 0(ma 7 ). Additional contributions in the two- 
photon-annihilation channel involve terms having a one-loop 
SE, vertex, or ladder correction on each side of the diagram, 
and two-loop self-energy, vertex, etc., corrections on a single 
side of the diagram. Calculation of these terms is in progress. 
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